Spin and Charge Texture around In-Plane Charge Centers in the CUO2 planes 
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Recent experiments on La2Cui_xLix04 show that although the doped holes remain localized near 
the substitutional Li impurities, magnetic order is rapidly suppressed. An examination of the spin 
texture around a bound hole in a Cu02 plane shows that the formation of a skyrmion is favored 
in a wide range of parameters, as was previously proposed in the context of Sr doping. The spin 
texture may be observable by elastic diffuse neutron scattering, and may also have a considerable 
effect on NMR lineshapes. 



< 
<N 

00 

o 

00 

o 



I 

o 
o 



Although it has been known for some time that only a 
very small concentration of holes doped into the antifer- 
romagnetic (AF) state of the cuprates suffices to suppress 
long range AF order, the close proximity to the metal in- 
sulator transition makes it difficult to distinguish clearly 
between the length scales of charge and magnetic distur- 
bances. Recently, experiments on the La2Cui_xLix04 
system where Li^ substitutes for Cu^"'" in the Cu02- 
planes show a similar rapid suppression of long range 
AF order but without a transition to a conducting state. 
[1] In this system clearly the hole charge introduced by 
the Li^-substitution is well localized, whereas the mag- 
netic perturbation must be long ranged. Several authors 
have in the past proposed that bound holes introduce 
skyrmion textures into the AF ordered state. [2-4] In 
this letter, we model the Li-substitution and examine 
this possibility in greater detail in the limit of dilute Li- 
concentration. Our conclusion is that skyrmion textures 
which are known to destroy long range magnetic order [5] 
accompany holes bound to charge centers, and we pro- 
pose that the skyrmion textures may be observable in 
neutron scattering and NMR experiments. 

The substitution of Li"*" for Cu^"*" in a cuprate magnetic 
insulator such as La2Cu04 has two effects. First, a Li"'" 
ion is a net negatively charged center, and secondly each 
Li"'" ion is accompanied by a doped hole. This hole in turn 
is attracted to the Li"'" site by a Coulomb potential, but it 
is not allowed to hop onto the impurity site. Additionally, 
there is a change of the singlet Cu^"'" state of the hole on 
neighboring Cu-sites due to a modification of the pd- 
hybridization on these sites. 

To gain an intuitive understanding of the effect of 
a partially localized hole on the local magnetic mo- 
ment at neighboring Cu^"'" sites in the plane, let us 
first consider the atomic limit. To zeroth order, a sta- 
tic hole which is attracted to a Li"'" impurity (located 
at site 0) by a nearest neighbor Coulomb interaction 
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where Gx and Gv are unit vec- 



sites) will reside on one of the Cu02 squares in the in- 
nermost "shells" - consisting of the 4 Cu02 adjacent to 
the impurity ion (Fig. 1). The energy of the resulting 
singlet is given by 
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where Eb is the binding energy of the singlet, while 
the factor a takes into account the anisotropy of 
the singlet due to its proximity to the impurity. 
Specifically, the distortion of the singlet structure 
((;/); = -i=(d|p^ — djpj^)) [7] around Cu^"'" sites adjacent 
to the impurity leads to a distorted singlet of the form 
$ = \[a(f)i -I- -^Vl - a'^{-(f)'2 + (f)?,- '/'4)], where the in- 
dex 1 refers to the 0-site adjacent to the impurity. 
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FIG. 1. Schematic plot of the 13-site cluster consisting out 
of 3 Cu-shells centered around a Li impurity. O-sites on the 
Cu-Cu and the Cu-Li links are not shown. 



tors linking the impurity site with neighboring oxygen 
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Introducing a hopping term, we find witliin second 
order perturbation tlieory about tlie atomic limit tliat 
tlie average liole occupancy on a Cu-site in tlie i"^ sliell 
around tlie impurity site is given by 



<n 



1 

>=4 



<5i,i(Ei 



Eb)2 + a,t2 



(El -Eb) 



5t^ 



(2) 



where t is the nearest neighbor hopping integral between 
the Cu-sites, and ai = 0, a2 = 4, as = 1. The correspond- 
ing magnetic moment in the i"^ shell is then given by 
< Si >= ^(1— < nf' >). [8] Thus, for a strong enough 
Coulomb attraction, the hole will be constrained almost 
entirely to the innermost Cu-shell around the impurity 
yielding an average magnetic moment of < Si >~ |, i.e. 
a reduction of 25%. In addition, corrections due to zero- 
point fluctuations are expected when exchange processes 
among the spins are included. 

Let us now turn to the question of the spin background 
created by a dilute concentration of holes, which are par- 
tially localized around impurity ions, in an otherwise AF 
background. While the motion of a hole about an impu- 
rity ion has strong quantum character, tending to spin- 
polarize the region around the impurity, the area far from 
the impurity may be treated in a semi-classical manner. 
[2] It is then natural to assume that the spin background 
will compromise between the (partially) polarized area 
around the impurity ion and the AF long-range order 
which is expected at infinite distance from the impurity 
by favoring a staggered skyrmion configuration. [2-4] In 
this sense, the "quantum" region where the hole moves 
rapidly around the impurity imposes a boundary condi- 
tion on the quasi-classical spin background, which in turn 
may be described by an 0(3) continuum model that is 
renormalized by quantum fiuctuations. [4] 

As first noted by Belavin and Polyakov, [5] skyrmions 
are a metastable solution of the 0(3) non-linear cr-model 
with Hamiltonian = J d^r(Vn)^. Their solution 
for the two-dimensional isotropic ferromagnet can be ex- 
tended to the AF case, [2,3] which - for a given choice of 
sign of the topological charge - has the form 
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where A is the core size of the skyrmion, and the alter- 
nating sign corresponds to the two sublattices. We will 
demonstrate that such a spin texture is indeed stable in 
the presence of holes that are partially localized around 
defects in the plane. 

In the following, the region around the impurity is 
treated exactly, i.e. by exactly diagonalizing a cluster 
of the kind shown in Fig. 1, while the semi-classical 
spin background enters the finite system calculation as 



an outer boundary condition through an exchange field, 
JS(r) = |n(r. A). Then the free energy of the whole sys- 
tem (finite cluster -|- surrounding field) is minimized with 
respect to A to deduce the optimal core size Aopt of the 



skyrmion field. A minimum at A 
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would imply that the system favors AF alignment, and 
spin currents introduced by defects would be irrelevant. 
On the other hand, a finite Xgpt corresponds to a stable 
skyrmion spin texture. 

We study the planar t-J model with one hole on the 13- 
site cluster shown in Fig. 1, with boundary conditions 
imposed by equation (3). The in-plane impurity site, 
located at the origin of the cluster, attracts the hole 
via a nearest-neighbor Coulomb potential (|V| « 0.5t for 
Li). The corresponding Hamiltonian is given by 

H = -t J2 (2L5j-+h.c.) 

<1,J>#0,(7 

+J J2 (Si-Sj-^n,nj)-V J2 °o°±x,y' (4) 
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where the c-operators are hole operators acting on non- 
doubly occupied states, J is the exchange integral, and 
t is the hopping amplitude. This system is not transla- 
tionally invariant because of the defect. Furthermore, the 
twist components (Sx,Sy) of the applied skyrmion field 
at the boundary of the cluster create non- vanishing ma- 
trix elements between different Sj^j subspaces of the en- 
closed system. In the absence of the outer spin field, the 
ground state of the impurity-containing cluster is dou- 
bly degenerate (S^^j = ±i), and has a non- vanishing as- 
sociated current, j« = it[(5L7-^a'Ci+i,<,' - ^l+i„T'a,a'C,,a'], 
where a=x,y,z labels the components of the current, and 
1 is a unit vector in the direction of the (anti-)clockwise 
motion of the hole around the impurity. [3] 

In the presence of the boundary condition imposed by 
the skyrmionic configuration of the far field, Sj^j is no 
longer a good quantum number. As shown in Fig. 2(a) 
for the parameter values J/t=0.4, V/t=0.5, and t=1.0, 
SjQj is a function of A which changes its sign around 
Xopt ~ 2 lattice spacings. In the texture form (Eq. (3), 
the sign of the sublattice magnetization is fixed at large 
distances (r A), and the z-component of the exchange 
field at the boundary of the inner cluster is either paral- 
lel (A ^R) or antiparallel (A ;^R) to this magnetization, 
leading to a net S^^j = ±i in the two limiting cases, A 
and A co, as seen in Fig. 2(a). 

A transition is also observed in the z-component of the 
associated current defined above. In Fig. 2(b) the cur- 
rent along the 8 nearest neighbor sites of the Li"'" ion 
is shown as a function of the far field boundary con- 
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At Xopt, it jumps from «-|-0.35t to 



-0.35t. 



Correspondingly, the total energy of the inner cluster. 
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shown in Fig. 2(c), has a minimum at A 



opt • 



Note that Eciuster('^ — ?> 0) = Eciuster('^ co) because of 
time-reversal symmetry. While the spin contribution to 
the total energy of the cluster is largest around Xgpt , the 
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magnitude of the Coulomb energy - which is about an 
order of magnitude smaller than the spin contribution - 
is maximal at A = and A = cxd, i.e. when the mobility 
of the hole is suppressed by the far field, and it is forced 
into the vicinity of the impurity. 

Special care has to be taken when including renormal- 
ization effects due to quantum fiuctuations in the spin 
background. [4] In particular, the « 40% reduction of 
the sublattice magnetization leads to a reduction of the 
classical spin stiffness = J • roughly by a factor of 
2. [10] A corresponding reduction also enters the applied 
spin field at the boundaries of the inner cluster. 

The energy of the skyrmion far field is given by 
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where ro(~ -^5) is the equivalent "radius" of the finite 
cluster which is treated exactly. Efar is thus a monoton- 
ically increasing function of A. The total energy of the 
system , including the inner cluster and the far field, is 
shown in Fig. 2(d). The effect of Efar is a slight reduction 
of the optimal skyrmion core radius Xopt- 

From Fig. 2, it is observed that the system favors a 
skyrmion texture with a core size of 2-3 lattice spacings. 
The results shown in this graph are found to be robust 
for a wide range of parameters (J,t,V), [11] leading us 
to the conclusion that a skyrmionic spin configuration is 
favored for the lightly Li-doped system. 
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FIG. 2. (a) Total magnetization, S^ot, as a function of the 
skyrmion core size A. (b) z-component of the current in the 
two innermost shells around the Li"*" ion. (c) Energy (in units 
of t) of the inner cluster (Fig. 1), including hopping, spin 
and Coulomb contributions, (d) Total energy, including the 
contribution from the far field (Eq. (5)). 

Possible experimental probes of a skyrmion ground 
state should examine its magnetic nature. The Fourier- 
transform of the staggered magnetization and the mag- 
netic structure factor exhibit a less sharply peaked behav- 



ior in k-space than ordered antiferromagnets. To quan- 
tify this statement, here we calculate these quantities for 
a small concentration of impurities, neglecting skyrmion- 
skyrmion interactions. [12] 

The Fourier-transform of the staggered magnetiza- 
tion, S? = (-ir4^ = (-ir4[l-j^A^], con- 
sists of a (5-function at Q = (tt, tt) corresponding to 
the AF order, and a tail due to the topological tex- 
ture, l[,5q,Q - :i^Ko(A|q- Ql)], where Ko(x) is the 
O"^ order modified Bessel function of the second kind. 
Hence, for a dilute concentration of skyrmions (with 
concentration Fimp), a deviation from the AF Fourier- 
transform of the staggered magnetization of the form 
rinip27rA^Ko(A|q — Q|) is expected. 

Similarily, elastic neutron scattering 

experiments, measuring the magnetic structure factor. 
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(SJSJ,), should also pick up a 



r,„,p4^A2Ko(A|q - Ql)] + 0(r? A4). (6) 



In Fig. 3, the magnetic structure factor is plotted for 
various values of A at a fixed impurity concentration Fimp 
= 0.2 %. 
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FIG. 3. Magnetic structure factor in the presence of a di- 
lute concentration of skyrmions with core size A, for an im- 
purity concentration Pimp = 0.2 %. A deviation from the 
AF signal ((5-function at Q = (tt, tt)) in the form of tails is 
observed. 

A skyrmion spin texture should also have an observable 
effect on NMR lineshapes. NMR measurements probe 
the response of the anisotropic hyperfine coupling to lo- 
cal Cu^"'" electronic spins and the isotropic transferred 
hyperfine coupling to the neighboring Cu^"'" spins. [13] 
If = AxSf X + AySf y + A.Sfz + B Sj is the hy- 

perfine field, then I(w) = ^^^^[^ — l^iP) gives the line- 
shape which can be observed in measurements of the zero 
field shift. In the presence of a small concentration of 
skyrmions this response function can be calculated from 
Eq. (3), assuming that at large distances from the im- 
purities the spins align antiferromagnetically within the 
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plane, i.e. z || la, where h is one of the in-plane axes. Re- 
sults for various small Li-concentrations using typical val- 
ues of the hyperfine couplings (Ax = Ay = 0, = — 4B) 
are shown in Fig. 4. While at concentrations < 0.06% 
the lineshape is observed to be rather sharply peaked at 
ujo, the resonance frequency for AF order, a shift of the 
peak towards higher frequencies, along with an asymmet- 
ric broadening, occurs at doping levels > 0.2%. From the 
corresponding continuum model it can be shown that the 
high energy cut-off for the zero field shift NMR signal is 
at 2wo- 



ity of the defect are reduced. A dilute concentration of 
Li-impurities is shown to favor a skyrmion spin texture, 
leading to "tails" in the static magnetic structure factor 
[12] and to a broad and asymmetric lineshape of the zero 
field NMR shift. 
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FIG. 4. Zero field NMR lineshape as a function of 
Li-concentration. The poles appearing in the response func- 
tion have been given a finite width of e = 0.05loo- 

From susceptibility measurements on impurity-doped 
cuprates it is known that the AF phase is destroyed by Li- 
substitution of approximately 3%, while Zn-substitution 
of order 30% is required to suppress AF. [1] Neither of 
these in-plane substitutions leads to a metallic ground 
state. In the case of a dilute concentration of Zn^"'" sub- 
stitutions for Cu^"'", no holes are introduced, and the spin 
background is only locally distorted by the presence of 
inert sites. Li"'" substitution, on the other hand, brings 
holes into the Cu02 planes, yielding a long ranged dis- 
tortion of the spin background (Eq. (3)). The holes 
are localized around the Li"'" ions due to attraction via 
Coulomb forces, since Li"'" has a net negative charge with 
respect to Cu^"'". [6] Out-of-plane Sr^"'" substitutions for 
La^"'" also attract holes via Coulomb interaction leading 
to a similar scenario for the spin texture. [2,3,14] How- 
ever, from transport measurements it is known that a 
Sr^"'" concentration of about 3% is sufficient to produce 
a metallic ground state. Thus we may conclude that in- 
plane holes are better localized by in-plane Li"'" than by 
out-of plane Sr^"'". 

In conclusion, we have studied the effect of a small 
number of in-plane Li"'" ions on the charge and spin struc- 
ture in Cu02 planes. Holes are found to be localized 
around these impurities, to which they are attracted via 
broken-bond minimization and Coulomb forces. Thus 
the local magnetic moments on Cu^"'"-sites in the vicin- 
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